Relaxation Phenomena in Ferrites 

By A. M. CLOGSTON 

(Manuscript received March 2, 1955) 

J. K. Gait has recently suggested a mechanism for the losses observed in 
ferromagnetic resonance and domain wall motion in single crystals of nickel 
ferrite containing small amounts of divalent iron substituted for divalent 
nickel. Each such substitution provides one electron able to move between the 
various octahedral sites in the crystal. Gait's suggestion is that the losses 
arise from a relaxation associated with the motion of these electrons. 

This paper develops a theory of this mechanism based on a thermodynamic 
model. Expressions are found for the velocity of domain wall motion, for 
the line width in ferromagnetic resonance and for the displacement of the 
field for resonance. 

INTRODUCTION 

In a recent paper, Gait ' 2 and Wijn and van der Heide have suggested 
a mechanism for the losses observed during ferromagnetic resonance and 
domain wall motion in single crystals of nickel ferrite having a small 
amount of divalent iron substituted for divalent nickel in the octahedral 
sites. It is supposed that each substitution of a divalent iron atom for 
a divalent nickel atom provides one 3d electron able to move from one to 
another of the various iron atoms on the octahedral sites. For conven- 
ience, we shall term such electrons the "free electrons". With the crystal 
magnetized in a particular direction, these free 3d electrons will always 
so arrange themselves as to minimize the free energy of the crystal. The 
free energy, therefore, is going to depend upon the direction of magneti- 
zation because of this mechanism, in addition to its usual dependence 
upon the crystalline and shape anisotropy. It is clear that there will be 
an additional torque acting on the magnetization. 

In what follows, we shall assume that the free electrons cannot instantly 
assume the equilibrium arrangements specified by the motion of the 



1 Gait, Yager and Merritt, Fhvs. Rev. 93, No. 5, p. 1119, 1954. 

2 J. K. Gait, B.S.T.J., 33, No. 5, p. 1023, 1954. 

3 H. P. J. Wijn and H. van der Heide. Rev. Mod. Phys. 80, p. 744, 1950. 
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magnetization but tend toward equilibrium with a relaxation time t. We 
shall show that this assumption results in energy being extracted from 
the motion of the magnetization and appearing as heat in the crystal. 
The above mechanism has been analyzed by Gait in reference 2 in 
connection with the velocity of domain wall motion in nickel-iron ferrite. 
A similar question has been analyzed by L. Neel 4 in discussing the motion 
of domain walls in iron. In Neel's case the relaxation is associated with 
the migration of carbon atoms between various sites in the crystal. 
The theory used here is essentially similar to that applied by N6el. 
We shall try, however, to emphasize more than does Neel the connec- 
tions with thermodynamics, and shall apply the theory not only to do- 
main wall motion but also to ferromagnetic resonance. The theory is not 
in accord with that proposed by Gait and the points of disagreement 
wall be mentioned below. 

1 . Kinetic Model 

We shall consider the assembly of N free electrons per unit volume to 
constitute a thermodynamic system in heat contact with the crystal 
lattice. We suppose that there are m possible sites per unit volume of 
crystal for the free electrons to occupy; m is, of course, just the number 
of iron atoms per unit volume lying in octahedral sites. We shall let 
£i(M) be the energy levels available to the electrons on the octahedral 
sites where M is the magnetization vector. We shall suppose that there 
are di levels per unit volume of energy £, and have therefore z2 d { = m. 
We shall let N> be the average number of electrons per unit volume in 
the level e< . If the magnetization is held steady in some direction, N> 
will approach an equilibrium value that will be denoted by N ira (M). We 
shall assume a universal relaxation time t such that AT,- approaches 
equilibrium according to the equation 

dNj = #,-. - Nj (1 x) 

dt t 

and this equation is to hold whether or not N ioo is a function of time. 
We observe that ^Ni = N and that therefore JZ dNt/dt = as must 
be the case. 

N i00 is obtained in the usual way by minimizing the free energy at 
constant temperature, and is given by 

tfn. = ZF^TTt (1-2) 

4 L. N(5el, Thdorie du Trainage Magndtique de Diffusion, Journal de Physique 
et Radium, 13, p. 249, 1952. 



N = E ^- £ //^_ L1 (1-3) 
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In what follows, we shall assume that the percentage of substituted iron 
is small enough so that rf,- ~2> N for each level. We have for the total 
number of electrons per unit volume 

di 

+ 1 

It is clear that e*'~ C/ must be much greater than unity for all i. Con- 
sequently, we may write 

N„ = N £^-=i^ (1-4) 

We now define the internal energy, the entropy, and the free energy 
per unit volume as 

U - E ft* (1-5) 

i 

S = k[N\nN -Eft In ft] (1-6) 

* 

F = U - TS (1-7) 

The equilibrium value of these quantities are obviously given by 

C/co = E # «•*< (1-8) 

5oo - k[N In JV - E ft- m #*J (1-9) 

Fco = Uao - TS* (1-10) 

If we take a time derivative of equation (1-5) we obtain 

dU v-* ., dEi . yr^ dNi /•. ,,\ 

*-?*<* +?«-* (M1) 

and in accordance with the first law, identify the rate of doing work on 
the system as 

while the rate at which heat flows into the system is 

dQ ^ dft , , 

2. Domain Wall Motion 

We shall first apply the kinetic model to a discussion of the velocity 
of a 180° domain wall moving perpendicularly to the (110) plane. In this 









742 



THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1955 



case the magnetization vector remains essentially in the (110) plane 
and its direction can be specified by the angle as in Fig. 1. We shall 
take the ?/-direction into the plane of the paper. Let us suppose that the 
wall is moving in the positive y direction, and that the magnetization 
turns through 180° from 0i to 2 as the wall sweeps by a given point. 



(001) 




(HO) 



Fig. 1 — Angles used in describing motion of domain wall in (110) plane. 

If we follow the procedure outlined by Kittel, 5 it is easy to show that 
the shape of a domain wall in the (110) plane for a crystal whose easy 
direction of magnetization is (111), is given by 



where 



A (|J = g(d) - g(eO 



g(0) = K\]4 sin 4 + sin 2 cos 2 0] 



(2-1) 



(2-2) 



K is the first order anisotropy constant (a negative number) and A is 
the usual exchange constant. We have sin 0i = \/%, cos 0i = ->/% so 
that 0(00 = K# and g(6) - ?(0i) - (-K/12)(2 - 3 sin 2 0) 2 . Equation 
(2-1) determines the wall except for its position along y. We have ig- 
nored in (2-2) a term depending on the magnetostriction which removes 
an ambiguity in the shape of the wall at the easy direction intermediate 
between 0i and 2 . This term is unimportant in what follows. 



6 C. Kittel, Physical Theory of Ferromagnetic Domains, Rev. Mod. Physics 
21, p. 656, Oct., 1949. 
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We shall now assume that the wall moves along y with a constant 
velocity v and maintains the static shape given by (2-1). At each point 
in the wall we can then write, 

dd dd tn 0> 

— = - v — (2-3) 

dt dy 

The domain wall is, of course, moving as the magnetization turns to 
line up with an applied external field H . The rate at which this field 
does work on the wall per unit area is given by 2H Mv. This quantity 
must be set equal to the net rate at which work is done on our thermo- 
dynamic system, or, what is the same thing, the net rate at which heat 
flows into the crystal lattice. We may write therefore, 

f°° dW 
2H Mv = ]_ ^-dy (2-4) 

In order to compute dW/dt from (1-12) we suppose that the domain 
wall is moving slowly and that N ix is therefore a slowly changing func- 
tion of time. Consequently we obtain from (1-1) an approximate ex- 
pression for Ni 

Ni -**• - T ~dT + T M + "--• (2 " 5) 

We obtain immediately from (1-12) to first order in t, 

The first term in this expression can be written as dF^/dt, and contributes 
to the integral over y a term 






dd 



We assume that F„ has cubic symmetry, and there is therefore no net 
contribution from this term. 

The second term may be treated as follows. Proceeding from (1-4) we 
find by differentiation 



dNi 



'00 



= -7^F Nu. ~ + 7^ ^ I E N+ § (2-7) 



dt kT ,0 ° dt kT ,M N f ,co dt 
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Then we have 

~ kT 2 ? N \_dt *t N dt] 

v dN<„ dEi _ _N ( \d h _ / & \T \ f *Y (2-9) 

V ~df It "" fcT * L^ d» J W 

where the angular brackets indicate a mean value over the equilibrium 
system. Upon substitution of (2-9) in (2-6) we find, 

dW dF m , N . TdSi , dEi .T . (dd\ 2 , m 

-dr = -dT + T kT { \_de- { W ) ] > W ( 10) 

By using the first law, we can write for the heat flow into the system 

dQ dU dF„ N .[da .(to.T./*V f2in 

It is clear that, in any process that takes the system between two iden- 
tical thermodynamic states, the terms dU/dt and dFJdt will integrate 
to zero. We must therefore regard the positive definite expression 

N , [d£i , dEi v T x 

T kf { l-de' { de ) ] } 

as an irreversible heat loss from the system, and it is just this heat loss 
that makes dW/dt greater than dF„/dt. 

Before proceeding we wish here to make certain comparisons with 
the theory used by Gait in discussing domain wall motion in Reference 2. 

8 It is possible to formulate this problem in terms of a temperature T K such 
that 



One then obtains 



where 



Ni = Ne- £ilkT "fL e- £ilkT * 



dt t 



Ot-'&M -<«» 
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From (15) to (24) in Gait's paper one finds, to first order in 



dW = T ft*« I dd 



dt dd 2 \dt 



which is to be compared with (2-10). This formulation seems somehow 
to neglect the equilibrium free energy of the free electrons, but since 
this term is indistinguishable from the ordinary crystalline anisotropy 
energy the omission is unimportant. The expression above, however 
must be interpreted as the irreversible heat loss and should therefore 
be a positive definite quantity as in (2-10). The quantity (d 2 gi ao /dd 2 ), 
however, is not positive definite, and this would seem to be a serious 
objection to Gait's formulation of the theory. In Reference 2, g ita is set 
proportional to g(6) defined in (2-2) so that {d 2 gi X /dd 2 ) in that case 
specifically has regions of positive and negative value, which should not 
be the case. 

Using equation (2-10) we find 



LJJ *- T 6£<[ar -< ^ > ]>(f)* 



2 N 

- fr BP 



(2-12) 



where (dd/dy) is taken as positive by convention. This integral is the 
net rate at which work is done on the wall and is positive definite. The 
corresponding expression developed by Gait (30) is 



£5*-*JC*®* 



We can only regard a positive value for this integral as accidental and 
to depend upon the shape of the domain wall. 

Combining (2-4), (2-12) and (2-1), we find for the velocity of domain 
wall motion 

2HoM 

V= (N\i/-*f* m t\** / d£ «\T\io 9o j i.ii, (2-13) 

This expression exhibits the characteristic dependence of v on 1/t dis- 
cussed by Gait in reference 2 and shown there to account for the very 
marked dependence of v on temperature found in his experiments. It 
also contains the obvious result that v depends inversely on the number 
of free electrons. 
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It would be very desirable now to proceed to calculate the coefficient of 
1/t. This can only be done by introducing some physical ideas to arrive 
at an explicit dependence of £,• on 0. This may be done according to two 
simple schemes that we shall consider later. First, however, we shall 
discuss the phenomenon of ferromagnetic resonance. One point that will 
soon become clear is that the present theory allows no direct comparison 
between domain wall motion experiments and ferromagnetic resonance 
experiments unless a particular model is adopted. 

3. Ferromagnetic Resonance 

We shall now apply the kinetic model to a calculation of the line 
width for ferromagnetic resonance and the corresponding displacement 
of the field for resonance. We suppose that the applied field and equi- 
librium direction of magnetization lie in the (110) plane at angles tj and 
do as shown in Fig. 2. The direction z lies along Mo . The direction x lies 
in the 1 10 plane perpendicular to z, while the y direction is perpendicular 
to the (110) plane in such a sense that x, y and z form a right-handed 
system. 

We begin by writing down an expression for the total free energy of 
the system. 

Ft = -ME cos (0 - n) cos* + 0(6, $) 

+ Z Nid - kT(N In N - E Ni In N t ) (3-1) 

i i 

where 6 and * are the instantaneous directions of the magnetization. 



(001) 



(010) 




(100) 



Fig. 2 — Coordinate systems used in discussing ferromagnetic resonance in the 
(110) plane. 
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It will be observed that and $ are not the usual polar coordinates but 
are angles especially suitable to describing events occurring in the (110) 
plane. The first term in this expression is the energy of the magnetization 
in the external field. The second term includes the crystalline anisotropy 
energy, the shape anisotropy of the specimen and any energy due to 
magnetostriction. The third term is the free energy of the free electrons. 
We now wish to obtain the x and y components of the torque acting on 
the magnetization. These torques will be given by 

dO d«P 

where the derivatives are taken with Ni held constant. This procedure 
is consistent with (1-12) for dW/dt. 
We obtain for T y , setting $ = 

T v = -ME sin (•-,)- 22 - £ N t ^ (3-3) 

do i do 

Returning to (1-1), we can write generally for JV, , ignoring a transient 
term, 

jV, = n - e-" T [ ' e" T ^2 dt (3-4) 

J dt 

If we insert (3-4) into (3-3) we obtain, 



T v = -MH sin (6-v)- -^ - Z N ia ^ 



(3-5) 



Let us expand this expression around the equilibrium angle 0o to obtain 

- (** " <* ~ «> + (§\ + ? [s (*<- £)]„) - (3 " 6) 



3eA /d-V,-.\ „ T /" ,,r de\ 

wL\rsrL § 1 e \M Ul 
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which can be expressed in terms of the equilibrium free energy as 

- (** o«(*-,) + (g), f +(§=)> (3-7) 

We choose O so that ME sin(0 o - v) + (dG/dd) Bo + (dFJdO), = 
and suppose that H is large enough so that Q is substantially equal to 
??. We shall now drop subscripts and suppose that the field and equi- 
librium magnetization are in the direction 0. We further recognize that 
A0 = M x /M. The expression for T v then becomes, 

+?(S( ! M'*/ v :« b « 

An exactly similar expression is obtained for T x 

-?(S)(30<~/ > *3* 

The torques T* and 5P y may now be substituted into the equation of 
motion of the magnetization 

™L= -yT-yMXh (3-10) 

at 

where 7 = ge/2mc and £ is the driving field. It will be observed that 
M x and M u will vary sinusoidally with frequency co of the driving field 
for small amplitude motion. In that case, we can write 

J dt 1 + (wt) 2 dt 1 4- («t) 2 

and similarly for M„ . It now becomes an elementary matter to compute 
the resonance line width and the field shift brought about by the re- 
laxation process. 



(3-9) 
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The total line width is given by, 

"-^?(ftX*) + ft)(*))n^ •* 

while the field required for resonance is increased by, 



\2 n i2 1 



8ff= -J-f^ + l^ 
2M \ dd 2 T a* 2 . 



2M V \\<W \ at? / T \a*/ \ a* // 1 



(3-13) 



+ (cor) 2 
The expression 



(ftX^ + ftXW) 



is very similar to an expression that arose in discussing domain wall 
motion. We can treat this term exactly as was done in arriving at (2-9). 
We shall find that 



(3-14) 



?((»)(*) + ft) (W) 

_N_ . VdEi _ . d£i v~] 2 , Yd£i _ . d£i .T. 
kT K lde ^ dd ; J ^ Ufa K d^ ) ] ) 

Let us first notice that AH and 8H depend upon the quantity 

which does not appear at all in (2-13) for domain wall velocity. Without 
making more explicit assumptions about e, , therefore, we cannot make 
any connection between the two experiments. The essential point here of 
course is that M is not constrained in the resonance experiment to move 
in the (110) plane as it is in the domain wall experiment. 

Let us now make the reasonable assumption that the quantity 

varies in the (110) plane with cubic symmetry. We write therefore that 

<Br-<y>I + BS-<S>I> (3 , 5) 

= V + QiH sin 4 6 + sin 2 6 cos 2 0] 
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When we come to discuss specific models, we shall give particular ex- 
pressions for p and q in terms of one unknown constant. We see that 
p > and p + q/S > to maintain the right hand side of (3-15) posi- 
tive. We observe that it is not possible to claim cubic symmetry for the 
terms 

®-<*>J> - d <[S-<3>]> 

separately. 

Using (3-15), the expressions for line width and line shift become, 

AH = F ffl [P + g(1/i Shl4 ° + Sln2 e C ° S2 6)] 1 +7cot)» (3_16) 

2MVa^ x a* 2 , (; , l7) 

- m [p + * ( " sin4 ' + sin2 * cos2 " )] rrfe 

Equation (3-16) for the line width shows a dependence on («r) that 
is characteristic of relaxation processes. If a measurement is made of 
AH versus temperature at a given frequency and in a given crystallo- 
graphic direction, a peak will be observed in the line width at very nearly 
the temperature where r = 1/w. The peak will not fall accurately at 
t = 1/w because of the temperature dependence of the coefficient of 
uvfy _j_ ( wr ) 2 but comes very close to l/o> because of the exponential 
dependence of r on temperature. Observations of this kind have been 
made by Gait, Yager and Merritt on nickel-iron ferrite and reported by 
them in Reference 1. 

If we may now extend the measurements over a range of frequencies, 
and observe how the line shape varies as a function of frequency at con- 
stant temperature, AH will have its maximum accurately at w = 1/t 
In principle, then, t can be determined as a function of temperature. 

If t is known at a particular temperature, a measurement of AH versus 
at constant temperature and frequency will serve to determine 
\/M N/kT p and l/M N/kT q at the temperature in question. 

Let us now consider the field shift bH given by (3-17). In the first 
place, this shift is not the ordinary displacement that accompanies the 
introduction of loss into any resonant system since such a field shift 
would be proportional to (AH) 2 . To the contrary, the field shift is of the 
same order of magnitude as the line width and we must consider seri- 
ously what effect this will have on the measurement of p-value. We have 
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assumed previously that F* has cubic symmetry, and in that case, the 
term in (3-17) depending on F x will be indistinguishable from the usual 
crystal anisotropy and need not be considered further. The second term 
we must consider more closely. 

The following expression is an identity. 

V + <lVA sin 4 + sin 2 cos 2 0] 

/ \ (3_18) 

= [V + |J " ^ (2 - H sin 2 - i% sin 2 20) 

Now, the second term on the right has the angular dependence in terms 
of which one usually expresses the first order anisotropy constant Ki . 
We may expect then that the measured anisotropy constant will have 
the form 

*' = *»-£ ©TO ««• 

where Ka is the static first order anisotropy constant. 

We may look, therefore, for a dependence of Ki on frequency and tem- 
perature like that given in (3-19) and make comparisons with the be- 
havior of the line width. 

If the total anisotropy field has been thus accounted for, we will be 
left with a line shift 

1 N ( q\ M* (3 _ 20) 

which will appear as a change in the ^r-value. If a measurement of the 
field H for resonance is made for a given co, T and 0, and if this field is 
corrected for the anisotropy field to give an effective field H e , we will 
have 

H -- y -mw( p + i)^tr' (3 " 21) 

We shall measure an effective 7, however, given by 

H e = - (3-22) 

Equating (3-21) and (3-22), we obtain 

or, if 7„ = g e e/2mc and Ag = g e — g, 
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*g = -J— JL'(p + £\ ( "?' (3-24) 

g 2H e M kT \ V T 5/ 1 + (cot) 2 

We, therefore, see that the ^-factor will also have a dispersion in fre- 
quency and temperature, and the measurement of ^-factor may indeed 
be the most sensitive way of observing the relaxation phenomenon. 

4. First Model 

There are two simple models we may adopt to obtain an explicit de- 
pendence of Ei on and $. In the first case, let us suppose that there are 
four distinguishable energy levels in the crystal £1 , £ 2 , £3 and £ 4 and 
that each is associated with one of the body diagonals of the cubic crystal 
according to the relation 

Ei = w cos 0» (4-1) 

where 0,- is the angle between the direction of magnetization and the 
i th body diagonal and w is some constant depending on temperature. 

This model is particularly tempting because there are four non- 
equivalent octahedral sites in the ferrite crystal, each associated with a 
given body diagonal. Expressed in terms of the direction cosines of the 
magnetization (0-10:203), we find 



Ei — 5- [1 + 2(ai<x 2 + 0:20:3 + 030:1)] 



£2 = o [1 + 2(aiQ:2 — 0203 — 030l)] 
O 



(4-2) 



£3 = 5- [1 + 2(— 010:2 + 0203 — a 3 ai)] 
o 



£4 = 5- [1 + 2(— aia-2 — 0203 + o 3 ai)] 
o 

while expressed in terms of the angles and <i>, the energies become, 
El = H! J 1 + [sin 2 cos 2 $ - sin 2 $ + 2\/2 sin cos cos 2 *]> 

f2 = H! ji + [ s in 2 cos 2 $ - sin 2 $ - 2\/2 sin cos cos 2 #]> 

f 3 = H! J 1 + [-sin 2 cos 2 * + sin $ + 2\/2 cos sin $ cos 4>] > 

= !? J 1 -f [-sin 2 cos 2 * + sin 2 $ - 2\/2 cos sin $ cos $] | 



(4-3) 



£4 = w 



RELAXATION PHENOMENA IN FERRITES 753 

Let us consider now how we shall calculate the various mean values 
that will be required. If we consider a set of values <7, we have 



(4-4) 



Z d ie - { «-w kT 

i 

We shall suppose that T is large enough so that (e ,• — (f, » <JC kT, and 
assume that rf, is the same for each level. In that case, approximately 

(9i)=lZg n (4-5) 

k „=i 

if there are k levels. We shall show shortly that this assumption is con- 
sistent with the experimental data. 

Using (4-5) we can calculate for our first model the following mean 
values, 

<B'-<£>]>=© ![1+W2 * (4 - 7) 

eft 



{ ldd " { de ; J + U* " ' ■'■'• ' 



= (l J [1 + 3 cos 2 20 + 4 cos 2 0] 



(4-9) 



= %w 2 [l - 2(K sin 4 + sin 2 6 cos 2 0)] (4-10) 

Referring back to (2-13), the velocity of domain wall motion becomes 

2HoM 



{ w)]/^C { ] :i : ' r, " : '- > " 1 --' A ^ n '"' 



v = 

■ ( 
.kT) V 12A Je, V3, 

2H<*M 



(4-11) 
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Comparing with (3-15) we find p = % w 2 and q = - l % w 2 . The line 
width for resonance becomes 

AH 

The anisotropy constant is given by 

*-*+»®« r rtfei (4 " 13) 

and the shift in {/-factor by 

Ag = 1 4 / tf \ , ^ (4 . 14) 

p # C M 15 \/c77 1 + (cot) 2 

We thus have four different measurements which we may try to fit with 
the same value of N/kT w 2 . It should be noted that (4-12) predicts a 
maximum line width in the (100) direction and a minimum in the (111) 
direction. 

5. Second Model 

For the second model, we shall assume that there are three distinguish- 
able energy levels d e 2 and e 3 each associated with one of the cube edges 
of the crystal according to the relation 

£i = w cos 2 di (5-1) 

where 0, is the angle between the direction of magnetization^ and the ith 
cube edge. Expressed in terms of the direction cosines of M we have 

fl = w «i 

£ 2 = w cm (5-2) 

&i = w a 3 
while expressed in terms of the angles and $ the energies become, 

£i = 77 (sin cos $ — sin *) 2 



€ 2 = ]j (sill cos $ + sin *) : 
e 3 = w (cos cos <£>) 



(5-3) 
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This model we now consider is the same as the one used by N6el in 
Reference 3 to discuss the relaxation of interstitial atoms in iron. Pro- 
ceeding as before we find, 

<s>-Gf) ^ 

<[§- <^>]^=2™ 2 sin 2 0cos 2 (5-5) 

<fe-<^>T>=^ 2 sin 2 * (5-G) 



(5-7) 



|_ </* d$ ' J 3 






3 
In this case we obtain, 

2H a M 



2 / 1 _• 4 n l -2 n 2 

= - w 1 - sin + sm cos 



v = 



(tL) A/t^ ( 2 ™ 2 sin 8 cos 2 I 2 - 3 sin 2 I <70 
■ * '' '^ (5-8) 



2H M 



A /EK(E 



w~ (.583) 



and by comparison with (3-15) p = and q = % w . The line width for 
resonance is 

*" - £ i © » e sw s + * • c ° s ' •) to (5 - 9) 

The anisotropy constant is 

The effective field at resonance is 



co 4 /AW 2 \ (cor) 2 

^ e " 7 " 15 VjfWy 1 + («r)» (5_11) 



and the {/-factor shift is 



v, i ^ (N\jMr (M2) 



<7 ZM/ \15/ \A-77 1 + («t) ! 
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Nio.75Fe2.26O4 after Yager, Gait and Merritt. 

In this case the line width due to relaxation is zero in the (100) direction 
and a maximum in the (111) direction. 

6. Comparison with Experiment 

We wish now to make some comparisons between the theory developed 
in the preceding sections and measurements upon nickel-iron fer- 
rite of composition Nio.75Fe2.26O4 made by Yager, Gait and Merritt. 
The data we shall consider were made upon a .009" sphere and are pre- 
sented in Fig. 4 and Table II of Reference 7. These data include a 
measurement of line width versus temperature for the (111) direction 
and (100) direction in the crystal; a measurement of the first order 
anisotropy constant (as Ki/M) versus temperature; and a measurement 
of the effective field at resonance, H e , versus temperature. It is observed 
from these data that the line width is maximum in the (111) direction 
and we shall therefore try to fit the data with the second model as dis- 
cussed in Section 5. This model requires the line width to be zero in the 
(100) direction. In fitting the data, we shall assume that the line width 
in the (100) direction represents non-angle-dependent contributions from 
other mechanisms and shall subtract this width at each temperature 
from the width in the (111) direction. The resulting curve is shown in 

7 Yager, Gait and Merritt, Phys. Rev., to be published. 
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Fig. 3. The data for Ki/M and H e are presented in Table I. The effective 
fields H e have been corrected in each case to a common frequency 24,388 
mc/sec. 

The equations we need are contained in Section 5 and are repeated 
here for convenience, 
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To proceed, we must make some assumptions about how w depends 
upon temperature, and shall therefore suppose that (Nw 2 /MkT) is a 
constant independent of temperature. If we assumed instead, for in- 
stance, that w was a constant, the curves to be calculated would differ 
considerably at low temperatures but not enough to affect the conclu- 
sions to be drawn. 

Proceeding on this assumption, we conclude from Fig. 3 that cor = 1 
at T = 145° and calculate that (Nvf/MkT) = 690. We may now use Fig. 
3 and (6-1) to calculate cor at each temperature. With this we may cal- 
culate H „ at each temperature, and the resulting curve is shown in Fig. 
4. Since we have no independent value for y, the curve has been fitted 
to the data at 160°K. It is seen that the predicted changes in H e are 
in the correct direction and are of approximately the right magnitude. 
The experimental points are not fitted in detail, however. The fit above 
85°K could be somewhat improved by assuming a different dependence 
of w on T. The theory as presently developed cannot explain the in- 

Table I — Anisotropy Field —Ki/M and Effective Field 

H,, for Various Values of Absolute Temperature 

After Yager, Galt and Merritt. 
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Fig. 4 — Measured and calculated values of effective field //„ as a function of 
absolute temperature. 
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Fig. 5 — Measured and calculated values of anisotropy field —K\/M as a func- 
tion of absolute temperature. 
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crease of H e at 4.2°K. Yager, Gait, and Merritt believe that the changes 
observed between 85°K and 4.2°K must be ascribed to some other 
mechanism. 

In arriving at (6-1) to (6-3) we have assumed as in Section 4 that 
w/kT « 1. At a temperature of 150°K, kT equal 2.07 X 10" 14 ergs. 
The molecular weight M of nickel ferrite is 234.39 and the density p is 
5.37 gms/cc. There are then pL/M = 1.38 X 10 22 molecules per cubic 
cm, where L = 6.025 X 10 23 is Avogadro's number. In nickel ferrite of 
composition Nio.75Fe2.25O4 one divalent iron is introduced for each four 
molecules, and the number of free electrons N is therefore 3.45 X 10" 
per cc. The magnetization of this ferrite is about 340 cgs units. We have 
determined above that (Nw 2 /MkT) = 690 and can therefore calculate 
that w/kT = 0.057 at T = 150°K. We have thus proceeded consistently 
in assuming w <3C kT. If w is constant, our calculation would not how- 
ever be correct below about 50°K or if w /kT is constant, below about 
10°K. It would be entirely possible to make accurate calculations for 
the entire temperature range and, by assuming special variations of 
w with temperature, perhaps account for some of the rise in H e at 
4.2°K. 

We must also recognize that the ratio m/N of levels to number of 
free electrons in this composition is 5. Since we have assumed the levels 
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Fig. 6 — The theoretical relation between anisotropy field —K\/M and effec- 
tive field H e compared with experimental points. 
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divided into three equal groups, di = %N so that we are hardly justi- 
fied in considering the electrons to be non-degenerate. We should then 
really use (1-2) instead of (1-4) in computing our averages. Further- 
more, when degeneracy effects are important, the relaxation equation 
(1-1) should be replaced by a more sophisticated expression. Since the 
general effect of degeneracy is to reduce the number of electrons taking 
part in the relaxation process, particularly at low temperatures, a calcu- 
lation considering degeneracy has a good chance of accounting for the 
rise in H t at liquid helium temperatures. 

In a ferrite of composition similar to that considered here, Bozorth, 
Cetlin, Gait, Merritt and Yager 8 report measurements of the static 
anisotropy which show Ki/M to be nearly constant between 300°K and 
196°K. If we make the assumption that the static anisotropy varies 
very little down to 85°K, we may use (6-3) to draw the curve in Fig. 5. 
It will be observed that the predicted changes in Ki/M are again in the 
right direction and of approximately the right size although the detailed 
fit is poor. 

We may make one further comparison with the data. By eliminating 
between (6-2) and (6-3) one obtains 

With the same assumption about K , a plot of {—Ki/M) versus H e 
should be a straight line with slope — }>£• The data has been plotted in 
Fig. 6 and a straight line of the appropriate slope drawn through the 
point at 160°K. Although only three experimental points are available, 
the fit is very satisfactory. 

The author would like to acknowledge many helpful discussions with 
J. K. Gait, P. W. Anderson, H. Suhl and L. R. Walker. 

8 Bozorth, Cetlin, Gait, Merritt and Yager, Physical Review, to be published. 



